Precise determination of the rfc mass and width in the 
radiative J/xj/ ^ r]cY decay 



Nora Brambilla*, Pablo Roig" and Antonio Vairo* 

* Physik-Department, Technische Universitat Milnchen, James-Franck-Str. 1, 85748 Garching, Germany 
^ Laboratoire de Physique Theorique, CNRS/Univ. Paris-Sud 11 (UMR 8627), Bat 210, Faculte des Sciences, 

F-91405 Orsay Cedex 

Abstract. We present an effective field theory based extraction of the rjc mass and width from a recent measurement by 
CLEG of the photon line shape in the 7/ 1/^ ^> rjcY decay. 

Keywords: Radiative decays of charmonia, quarkonia, non-relativistic effective-field theories, QCD 
PACS: 13.25.Gv,14.40.Pq,12.38.-t,11.15.Me 



We study the radiative decay J/\j/ —i' rjcY within an 
effective field theory (EFT) framework, namely poten- 
tial non-relativistic QCD (pNRQCD) [1, 2, 3]. Our mo- 
tivation is to take advantage of the new measurements 
made by the CLEO [4] and KEDR [5] collaborations 
of the photon line shape, in order to obtain precise de- 
terminations of the rjc mass and width (for a recent re- 
view see [6]). Magnetic dipole transitions in quarkonium 
were studied in pNRQCD in Ref. [7], where, for the 
— > rjcY decay, a branching ratio consistent, within 
errors, with the PDG value [8] was found. We apply the 
same formalism here. Before discussing it, however, we 
briefly review the CLEO and KEDR analyses. 

To fit its data [4], CLEO uses two background sources, 

• a Monte Carlo modeled background for spurious 
J/\lf with shape 



bkgi (By) = N f e-5.720£, ^ ^Q_44^ g-33.567£, 



and a freely fit background for J/xj/ ^ ^'^X and 
non-signal J/y/ ^ X with shape 

bkg2(£y) =A+BEy + CE^, 



and a theoretical Une shape given by 

flieory(£'y) — x BWrei(£'y) x damping(£'y), 
where the relativistic Breit-Wigner distribution is 



BWrei(£y)-' = (Mj/^-2Mj/^Ey-Ml 



2M,/,„Ey] r 



Mji^ and M-^^ stand for the masses of the y/y/ and Tj^ 
respectively, Y-q^ for the rj^ width, and 



damping(£'y) — e 



-£2/(8/32) 



The damping function accounts for the overlap of the two 
quarkonium states, assumed to be described by wave- 
functions of an harmonic oscillator, and the photon. The 
natural scale of j3 is then the typical momentum transfer 
inside the charmonium ground state, which is about 700 
to 800 MeV. This implies almost no damping for photon 
energies smaller than 500 MeV, which is consistent with 
the multipole expansion of the electromagnetic fields. 
The value that comes from CLEO's fit is, however, an or- 
der of magnitude smaller, /3 = (65.0 ±2.5) MeV, which 
implies the vanishing of the signal for photon energies 
of few hundred MeV. The CLEO analysis yields the val- 
ues M^^, = (2982.2 ±0.6) MeV and r,,^, = (3L5±L5) 
MeV. 

The Tjc line shape is also studied by the KEDR collab- 
oration in [5]. Their analysis is similar to CLEO's one, 
but an alternative damping function is used as well; 



damping' (iiy) = 



p2 

^peak 



EyEpQ^ 



(Ey -Epeak) 



where fipeak is the most probable transition energy. If 
CLEO's data are used, the analysis of KEDR gives 
Mrf^ = (2982.4±0.7)MeVandr,,^, = (32.5 ±L8) MeV 
when fitting with damping(£'y), and M-q^ = (298 L8± 
0.5) MeV and T^^ = (33.6 ± L9) MeV when fitting 
with damping'(£'y). However, an analysis of KEDR's 
own preliminary data gives different values: M-q^. = 
(2979.7 ±L6) MeV and = (26.9 ±4.8) MeV when 
fitting with damping(£'y), and M^,, = (2979 .4 ±L5) 
MeV and L^,. = (27.8 ±5.1) MeV when fitting with 
damping' (Ey).' The discrepancy between different val- 



' Incidentally, the KEDR results are closer to previous rj^. mass mea- 
surements from and ^{2S) — > T/r7 decays (averaging (2977.3 ± 
1.3) MeV [9]), while the CLEO value is in agreement with the re- 



ues of the tJc mass and width is larger than the exper- 
imental sensitivity, which highlights the importance of 
performing a critical analysis of the theory inputs. 

Potential NRQCD exploits the hierarchy of scales in 
the problem and allows to express physical observables 
as systematic expansions in the ratio of these scales. 
Heavy quark-antiquark bound states are characterized by 
a number of scales: the heavy-quark mass, the typical 
momentum, (/?), exchanged by the quarks (this is also of 
the order of the inverse of the typical size of the bound 
state, !/(?")), the binding energy, the typical hadronic 
scale Aqcd, and possibly other smaller scales. In the 
transition J/Xj/ ^ ricY, the relevant scales are the charm 
mass OTc. which is much larger than the next-to-largest 
scale, which is (p) ~ 1 / (r) ~ 800 MeV, which in turn 
is larger than Aqcd- The binding energy of the J/\lf is 
Ej/y, ~ 500 MeV, which is smaller than 1/ (r). There 
are also two smaller scales that have to be considered 
in addition: the hyperfine splitting, which is M//^ — 
Mf)^ ~ 120 MeV and smaller than Ej^^,, and the width 
of the Tjc, Fjj^ 30 MeV, which is smaller than the 
hyperfine splitting. In the radiative decay J/yf Xy 
around the J]^ peak, we consider photon energies that 
vary between MeV and 500 MeV< 1/ (r). Under these 
conditions, we may describe the charmonium ground 
state in weakly coupled pNRQCD (because the system 
is non-relativistic and the typical momentum transfer 
is larger than Aqcd). couple the electromagnetic fields 
to pNRQCD and multipole expand the electromagnetic 
fields (because the photon energy is smaller than the 
typical momentum transfer in the bound state). 

Three main processes contribute to J/^ ^ Xy within 
pNRQCD in the energy range of interest (0 MeV< By < 
500 MeV): 

• the Ml transition J/xj/ ^ rjcY ~> Xy, 

• the El transitions J/y -)• Ze(o,2) (l^)7->- -^7. 

• fragmentation and other background processes in- 
cluded in the background functions. 

The evaluation of the magnetic dipole contribution 
yields (see Fig. 1) 



dE. 



(1) 

It has been pointed out in Ref. [4] that the dependence on 
Ey is responsible for the asymmetric shape of the photon 
spectrum. 



suits obtained in yyfusionand production (averaging (2982.6±1.0) 
MeV [9]). 
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FIGURE 1. Ml contribution to the 77c line shape. The black 
dot stands for four-fermion operators that contribute to the 
decay width of the intermediate state. 




FIGURE 2. El contribution to the tj^ line shape. 



The electric dipole contribution is (see Fig. 2) 

2 

448 - ^■'/'^^^^ 



dr^i 



dEv 



-a 
243 OTc 
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-a. 



mi 



ae(Ey) 



(2) 



where 0//i(/(O) is the 7/ wavefunction at the origin. The 
function a^{Ey) has been discussed in Refs. [10, 11]; a 
closed analytical form has been derived in [12J. 

In the weak-coupling regime, the typical momentum 
transfer in the charmonium is of order nicO^, the bind- 
ing energy of the ground state is of order nicd^ and 
the hyperfine splitting is of order nicaf. The magnetic 
and electric dipole contributions are of equal order for 
lUc Os ^ Ey^ mcOC^. The magnetic contribution com- 
pletely dominates the electric one in the peak region 
(rric (X^ ^ Ey ^ rric a^) and it also dominates by a fac- 



tor 



-J|^,I (Mj/v-^Vc) ~ ^/(^s for Ey < nic at 
We use as theoretical Une shape the sum of Eqs. (1) 
and (2), and as background the sum of bkgi(£'y) and 
bkg2(£y). As predicted by the power counting and nu- 
merically confirmed, for photon energies < 500 MeV, the 
electric dipole contribution is negUgible with respect to 
the magnetic one. The signal shape has been convolved 
with a Gaussian resolution function, whose resolution 
width is 4.8 MeV [4]. Our best fit is shown in Fig. 3. The 
fitting parameters are M,)^, r,j^, an overall normalization, 
the signal normaUzation, and the background parameters 
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FIGURE 3. Our best fit to CLEO's data for the photon spec- 
trum in — s- Tjc7 using Eqs. (1) and (2) for the theoretical 
signal together with the two background sources bkgj (Ey) and 
bkg2(£y)- 

A, B and C. In particular, the best fit line shape parame- 
ters for the rjc are (errors are only statistical) 

Mn^. = (2985.9 ± 0.6) MeV, 
r^^, = (28.6 ±0.2) MeV. 

The main differences between our analysis and the one 
of CLEO are summarized in the following. 

• At leading order in the multipole expansion of the 
photon field, which is justified in the range Ey < 500 
MeV, and in the non-relativistic limit, the overlap 
integral in the magnetic dipole transition is equal 
to one. This amounts to setting equal to one the 
damping function in CLEO's analysis. We have 
checked that higher-order terms in the multipole 
expansion are indeed negligible. Our conclusion is 
that CLEO's damping function has no theoretical 
justification as an overlap integral, which is also 
signalled by the unnatural scale of the parameter 
j3 . The absence of a damping function accounts for 
about 50% of the difference in the rjc mass determi- 
nation in our analysis and in that one of CLEO. 

• Another difference is that in our analysis the dif- 
ferential width turns out to be proportional to a 
non-relativistic Breit-Wigner distribution (xEy), 
see Eq. (1), while CLEO uses a relativistic Breit- 
Wigner distribution. The relativistic Breit-Wigner 
distribution resums some classes of relativistic cor- 
rections, without including others of the same size. 
Partial resummations could lead to spurious effects, 
while the use of a non-relativistic Breit-Wigner dis- 
tribution appears to be more justified in the frame- 
work of a systematic relativistic expansion, like 
the one provided by pNRQCD. This difference ac- 
counts for about the remaining 50% difference in 
the determination of the Tj^ mass. 



In summary, radiative decays of quarkonia and specifi- 
cally the transition 7/1/ —> TjcY may be investigated in an 
EFT framework that systematically exploits the hierar- 
chy of scales in the system (pNRQCD). The total transi- 
tion width has been calculated in [7] including the next- 
to-leading order relativistic corrections. Within a large 
theoretical uncertainty that determination is in agreement 
with the experimental value. In this work, we consider 
the photon line shape in the non-relativistic limit. We ob- 
tain a best fit, which is in good agreement with CLEO's 
experimental determination, see Fig. 3. However, theo- 
retical errors have not been included so far in our analysis 
and we expect them to be larger than those coming from 
the fitting accuracy. For instance, relativistic corrections 
could impact the rjc mass by coiTections as large as those 
induced by the difference between the relativistic and the 
non-relativistic Breit-Wigner distribution. Under investi- 
gation is also the extraction of the J/xj/ ^ Tj^. 7 branching 
ratio from the photon spectrum. 
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